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Abstract. We present a functional renormalization group approach to the 
zero bias transport properties of a quantum dot with two different orbitals and 
in presence of Hund's coupling. Tuning the energy separation of the orbital 
states, the quantum dot can be driven through a singlet-triplet transition. Our 
approach, based on the approach by Karrasch et al. |15| which we apply to 
spin-dependent interactions, recovers the key characteristics of the quantum dot 
transport properties with very little numerical effort. We present results on 
the conductance in the vicinity of the transition and compare our results both 
with previous numerical renormalization group results and with predictions of 
the perturbative renormalization group. 
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1. Introduction 

Quantum dot systems are ideal laboratories to study quantum-impurity models and 
allow to investigate many-body effects which result from the coupling of the quantum 
dot to a Fermi bath in the leads. If the ground state of the isolated quantum dot is 
degenerate, the coupling to the Fermi sea of the leads can induce a coupling of these 
different isolated ground states via virtual transitions and can yield a new strongly 
entangled many body state. A famous example is the 5=1/2 Kondo effect in which 
the leads completely screen the spin on the quantum dot and perfect transmission 
across the dot is achieved. pQ Kondo physics have been suggested to play an important 
role in a variety of the mesoscopic transport phenomena, in particular in formation 
of the famous "0.7 anomaly", where the quantized ballistic conductance through a 
quantum point contact acquires an additional conductance step at G — 0.7Go, Go 
being the quantum of conductance. [2] The 5 = 1/2 Kondo effect can arise, if an odd 
number of electrons are trapped on the dot. If an even number is trapped, a degenerate 
ground state can arise if a Hund's coupling leads to a finite spin of the ground state. If 
we consider a dot with two orbitals, separated by an energy difference S, and occupied 
by two electrons in the ground state, then a transition from a non-degenerate spin- 
singlet to a threefold degenerate spin-triplet ground state can be achieved by adjusting 
the energy difference of the two orbitals and/or the Hund's coupling strength. This 
singlet-triplet transition has been studied with a variety of techniques and can be 
detected directly by a peak of the finite temperature conductance at the transition 
point. OS] For a review on the transport properties of two-electron quantum dots, see 
e.g. 0. 

A robust and reliable technique for studying few level systems, such as quantum 
dots or nanotubes, is the numerical renormalization group (NRG) technique, [5J [7] 
which has also been applied in investigations of the singlet-triplet transition. [51 151 ITT] 
It is however limited to relatively simple geometries since NRG calculations for systems 
with more than two dots are numerically prohibitive. 

An alternative approach is the functional renormalization group (FRG). At the 
moment, available truncation schemes of the FRG equations allow for an accurate 
calculation of the dynamical properties of the single impurity Anderson model only 
at weak coupling, see Refs. [Jj3 [H]. In the limit of a vanishing bias voltage 
and at zero temperature, a calculation of the transport through the quantum dot 
requires however only knowledge of the local Green's function at the Fermi energy of 
the leads. A rather simple truncation scheme of the FRG was recently proposed in 
Ref. |15j . which, despite its simplicity, is accurate not just at weak coupling but also at 
moderate and even at rather strong coupling. The advantage of this approach, which 
we refer to as the static approximation since the frequency dependence of all vertex 
functions is neglected, is that it is considerably faster than NRG, and also simpler than 
other alternative approaches such as Fluctuation Exchange Approximation (FLEX) 
[T5] or field theoretical RG approaches [T7j and can be easily applied also to multidot 
geometries. [T5] The approach in Ref. [TS] was devoted to the study of density-density 
interactions and a Hund's exchange was not discussed. Here, we use the approach to 
investigate QD's with Hund's exchange. We present results for a two level QD with 
Hund's coupling and compare our results with previous works. 
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2. The model 

Considering a QD with two orbitals, there are multiple interactions to account for. 
Firstly, there are two intra- and one inter-orbital density-density interaction terms. 
Secondly, a Hund's spin interaction can exist between the two orbitals (in the static 
approximation, there is no independent intra-orbital spin interaction since for 5 = 1/2 
it can always be written as a density-density interaction). We thus arrive at the 
Hamiltonian of the isolated dot of the form 

Hd = ^2 \^A<j(^ Aa aAcr + E Ba a) Ba a B <j 

a 

+ Uab {p>A ~ - 1) + U A (nAt - \)( n Al - ^) 

+ U B (n m - ^){n B i - i) 

+ - ^2 Vo-fi ■ °V 3 o-4 a AtT 1 a Aa 2 a B a 3 a Bcr i , (1) 

0\ ,CT2 ,(T3 ,CT4 

where Ua and U B are intra-orbital interaction energies, Uab the inter-orbital 
interaction energy and J is the Hund's coupling energy. Here, cr = (cr x , a v , o~ z ) are 
the Pauli matrices with elements cr rTia2 . In the presence of a small magnetic field H 
the energy levels are Zeeman-split, 

E Aa =E A +ah, (2a) 
E Ba = E B +ah, (2b) 
with a = ±1, h = gfi B H/2 (fi B is the Bohr-magneton and g the electron g factor) 
and 

E A = V g -5/2, (2c) 
E B = V g + 5/2, (2d) 
where V g is the gate voltage and 8 is the difference of the two orbital energy levels. 

Note that V g — corresponds to half filling of the dot, as we put the Fermi energy of 

the leads to be zero. The leads are described by 

Hr + H L = ^ £k (c RT<jk CR Tak + c{ T ak C Lrakj , (3) 
rak 

where k labels the momenta in the leads and r = A, B denotes the different orbital 
quantum numbers. The coupling of the leads to the dot is described by the tunneling 
term 

HT = J2 (*«r C L-<7fe°™ + l ^ C \rak a r<y + kc) • (4) 
kra 

Note that the hopping conserves the orbital quantum number t — A, B, which makes 
the model appropriate for vertical coupled dots. In contrast, in lateral quantum dots, 
the orbital quantum number is not conserved. Lateral dots are described by a different 



tunneling term of the form Ht = Sfco-T y'R c Rcrk a To- + iZ/ c L CT fe a TCT + h.c.J and both 

orbitals share the bath electrons H R + H L = J2 a k £ k [ c Rak c R<^k + c \ak c L<yk*j ■ The 
absence of the orbital quantum number in the leads allows for virtual processes where 
the t quantum number of the dot electrons are changed. Later quantum dots where 
studied e.g. in Refs. [HI M EH HI]- Vertical quantum dots have been investigated by 
Pustilnik and Glazman |18| using scaling techniques and Izumida and Sakai [8 using 
NRG. 
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Conductance calculation 

Since the tunneling Hamiltonian Eq. Q conserves the orbital quantum number and 
the spin projection along the magnetic field, the total (zero bias) conductance through 
the dot is the sum of the partial conductances through each orbital level and spin 
projection. These can be obtained from the single-particle Green's function using the 
Landaucr-Biittikcr formula, [19 

G = G Ta , (5) 

with 

e 2 

G T a = J-Ptt^Lt^RtI Lt + T ' Rt ) ■ (6) 

Here, /9 Tff (0) is the density of states of the quantum dot at the Fermi energy and 
po is a broadening of the level with p being the density of 
states in the leads. The dot density of states can be obtained from the Fourier- 
transform of the retarded Green's function of the dot in presence of the leads, 
Gra(t-t') = -ie{t-t'){[a Ta {t),aX^)}), 

P r a {e) = --lmg Ta {e + iQ+). (7) 



3. Singlet-Triplet transition 



Here we concentrate on the situation with Ua = Ub = Uab so that, for total 
occupancy N = Yl T a nTcr ~ ^, there is no capacitative energy associated with 
transferring an electron from orbital A to orbital B. There is however a total energy 
cost which arises from both the Hund's coupling and from the energy difference of 
the orbitals. If both electrons occupy the lower A orbital, they necessarily form a 
(local) singlet. On the other hand, if one electron occupies the lower A orbital and 
another one the upper B orbital, then for J < the lowest energy can be achieved for 
a ferromagnetic spin alignment of the two electrons. The difference between the local 
singlet and the triplet state is 

K = S - J/4 , (8) 

and we thus expect a singlet-triplet transition at K = 0, where the configuration 
changes from a non-degenerate singlet to a threefold degenerate triplet state. In 
reality, because of renormalization effects, the transition is shifted from the point 
Kq = 0. One can however replace Kq by the renormalized quantity K such that 
transition takes place at K = 0. 



Symmetric coupling 

Here we focus on the symmetric situation where t^ = t^, tji T = tn and thus 
t 2 = t\ T + t 2 Rr is independent of t. On the triplet side of the transition the spins 
are locked together in a spin-1 state and the problem becomes a two channel spin-1 
Kondo problem with a completely screened groundstate and maximal conductance 



Functional renormalization group approach to the singlet-triplet transition in quantum dotsb 



at T = 0, while on the singlet side, conduction is effectively blocked. The analysis 
of Ref. |18j . based on Anderson's poor-man-scaling technique, predicts a monotonous 
decrease of G with K on the singlet side K > 0, 

G/G = Bin- 2 [K/T ] , B = ^-^j (10) 

with Go = 2e 2 /h,X = 2 + V5 and the characteristic energy scale 

T ~ f7exp[-r 7rLT/8r] (11) 



where To ~ 0.36. The asymptotic behavior (10) is valid for K 3> To. At finite 
temperature, the conductance develops a peak near K « 0. Both the asymptotic 
behavior and the peak near K = are recovered also in the FRG approach discussed 
below. 



4. Functional Renormalization Group 



According to Eqs. ( 6j7 ) the calculation of the conductance relies on the retarded 



Green's function, which we evaluate using the Matsubara imaginary time formalism. 
The full Green's function is defined as 

[Graiie)]' 1 = [GoM^V 1 - ^r*{ie). (12) 

where E r(T is the particle self-energy and Go,t<t is the noninteracting Green's function 
of the dot which, after integrating out the leads, has the form 

Go,ra(ie) = . — , (13) 

is + il T sgn(e) - h Ta 

where T T — 2irt 2 pv, t\ = |£lt| 2 + \tR T \ 2 and po is density of states in the leads which 
we take to be independent of r and a. All interaction effects are contained in the 
self-energy and below we calculate it using the non-perturbative FRG technique. 

The functional RG approach is based on an exact flow equation of the generating 
functional of irreducible vertices which can be used to derive flow equations for the 
irreducible vertices themselves. [2"U1 |2~T1 12"2"] The irreducible two point vertex is the self 
energy, the irreducible four point vertex corresponds to the fully renormalized two- 
particle interaction, and so on. This yields a hierarchy of coupled flow equations for all 
vertices which we truncate at the four point vertex, assuming that all vertices of higher 
order do not play an important role. In the scheme we use below, both the two-point 
and the four-point vertex are further assumed to be independent of energy, a rather 
drastic approximation which however was shown to yield very accurate results for the 
zero temperature linear conductance of coupled quantum dots. t 15i In the geometry of 
the quantum dots studied here, the two-point vertex is further diagonal in both orbital 
and spin quantum numbers if we choose the spin quantization axis to be parallel to 
the direction of the magnetic field. 

The flow equations emerge from introducing an infrared cutoff A in the non- 
interacting Green's function such that Gq — > for A — > oo and Gq = Go for A = 0. 
We introduce the cutoff via a multiplicative step-like regulator function 

Gt a (e) = 0(\e\^A)Go^, (14) 

where 9 is the Heaviside step function. Taking the derivative of the equation for 
the generating functional with regard to A then yields differential equations for the 
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irreducible vertices, and integrating from A — > oo down to A = yields the fully 
renormalized vertices. 

The flow equation of the self-energy £ A = £ A 1(Tl (the superscript A indicates that 
it is the self-energy at the cutoff value A, this applies to other quantities as well) takes 
the form 

9 ^ = E / £^ 0+ ^ A ( £ )7a 4) (1,2;2,1) (15) 

where 7 ^ 4 ' ) is the irreducible four-point vertex and the lower indices 1, 2 refer to both 
orbital and spin quantum number. The single-scale propagator is defined as 

G£(e) = -[Gi(£)] 2 M , 3oM~ 1 > ( 16 ) 

where G^(e) = £ A riffl (e), with £ A TlCTl (e) defined in Eq. 

Our sharp cutoff regulator yields a simple structure for both the single-scale 
propagator and the full Green's function, with 

qa ( p \ eq £ |-A) 

where the initial level energy and the self-energy have been merged 

Et = E Ta + £ A CT . (18) 

With the approximation of an energy-independent four-point vertex, its flow is 
determined by the equation (we neglect here the contribution of the six-point vertex) , 

<rf(i',2';2,i) = - / g£{e A m A (- £ ) 

J 3,4 

x 7 i 4) (l',2';3,4) 7 i 4) (4,3;2,l) 
+ 2S 3 A ( £ )S 4 A ( £ ) 

x [ 7 i 4) (l',3;4,2) 7 i 4) (2',4;3,l) 
- 7 i 4) (l',3;4,l) 7 | 4) (2',4;3,2)]}. (19) 

This flow equation involves products Q(x)S(x) which must be evaluated according to 
the prescription f[Q(x)]S(x) — > 5(x) f Q dyf[y] with some function f[y], see Ref. [21]. 

With the approximation of energy-independent vertices, the conductance 
contribution from each re, Eq. (J6|, takes a particularly simple form, 

e 2 r 2 

G ™ = Jfl^Ef„' (20) 
where we further assumed that = t^ T and where E ra = \\m E^ a 



Functional renormalization group approach to the singlet-triplet transition in quantum dots! 




fcCc( 



+ 2 





tf W 



; +24) +2+ 



.a 





+ 




Figure 1. Flow equation diagrams of the four-point vertices and the self- 
energies. Dashed and full lines are the propagators for the r = A and r = B 
orbitals, respectively. White (filled) 4-armed circles are the intra-orbital (inter- 
orbital) four-point vertices. Propagators with a vertical dash denote single-scale 
propagators defined via Eq. l |16| l and the circular dot on the l.h.s. denotes 
differentiation with respect to the cutoff A. 



4-1- Four-point vertices 

The frequency-independent intra-orbital interaction is completely specified by the 
interaction constants V a an d Ub- Furthermore, in absence of a magnetic field and 
of a direct hybridization of the two orbitals, the SU(2) spin symmetry allows only 
two types of inter-orbital interaction: a charge interaction which is specified by Uab 
and a Hund-type spin interaction of magnitude J. We can thus write the properly 
anti-symmetrized four-point vertex as follows, 



(4) 

7a 



„, inter 

7a 



„ Intra 

'7A : 



(21) 



with (S (Tl a2 denote Kronecker deltas) 



inter 

7a 



intra 

7a 



t=A,B 



.ttA 



(22a) 



(S Ti aS T2 B + S Ti bS T2 a) I — Y~ {^aia'^a 2 a'^TxT[K^ 



Functional renormalization group approach to the singlet-triplet transition in quantum dots8 



- a (Tl(T' 2 ■ £r 0- 2 CTj<W^T2T : ;) j ' (22b) 

where we dropped the arguments of the vertices, e.g. we implied ^ tor = 
7 A ntor (l', 2'; 2, 1) etc. These are all possible static two-particle correlations in presence 
of SU(2) symmetry and for two orbitals, in absence of a single particle hybridization 
of the two dots. 

4-2. Flow equations 

Fig.[l]shows a diagrammatic representation of the flow equations for the two-point and 
four-point vertices. The flow equation of E^ a directly follows from the flow equation of 



the self-energy Eq. ( 15 1, together with the form of the four-point vertex in Eqs. (22a 



and (226). We keep here a er-dcpendcncc of the self-energy which arises in presence of 



a magnetic field. We find 



9aEL = 



(E4)2 + (A + r r ) 2 

, E^(Ui B + J A /2)/(2w) 



(£7A)2 + (A + r,)a 

(23) 



(U£ B - J A /2)/(2n) 



The flow equations of the coupling constants which parametrize the four-point vertex 



follow from Eq. (19 1. For small magnetic fields, we neglect the effect of the magnetic 
field on the four-point vertex, since these yield corrections in the self-energy of only 
second or higher order in the magnetic field. We thus replace the cr-dependent energy 
levels £y CT by the average value Etf = {E^ + E^)/2 in the flow of the four-point 
vertex and arrive at 

5a£/ t a 



AttA\ 



2 (TO , 
7r[(£A)2 + (A + r r ) 2 ] 2 

?A\2 IK i -r \2l r3 



| mf-{K + T f f][{j^-{U^ B f 
2^[(^)2 + (A + r,)2] 2 

o rrA ^ [(E?) 2 -(A + Tr) 2 ]U T A Ui B 

okU ab = - 2^ r— — — - — — — 72 — 



(24a) 



t=A.B ^[(W + (A + r T )2] 

[f. 



[3E*E% - (A + r A )(A + r B )][f j a 2 + (u£ B f 



9\Ja 



*c[(Efi* + (A + r A )*] [(^ + (A + r B ) 
[(£ r A ) 2 -(A + r T ) 2 ]i/ A j A 

^— ■> _IVt^A\9 , IK , N2l 2 



2 



(246) 



~=a,b 7r[(iJA) 2 + (A + r r 



, 2E\eUaU^b - (A + r A )(A + TbW/2 

n[(E^ + (A + r A y}[(E^ + (A + T B y] ■ 1 > 



where r — B for r — A and vice versa. 



5. Results 



We have solved the flow equations ( 23|24c ) numerically, starting from a very large 
cutoff A and integrating the flow equations to A = 0. As discussed in Ref. [TS], the 
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appropriate initial values at Ao — > oo of the flow parameters are 

=E Tryi (25a) 

= U T , (256) 

U% = U AB , (25c) 

J Aa = J. (25 d) 

In particular, the bare one-particle contribution arising from the Uab and Ua/b terms 
in Eq. ([I]) are exactly canceled when integrating the flow equations from A = oo to a 
large but finite Ao.[TS] We limit the discussion to cases with a ferromagnetic coupling 
J < 0, for which the singlet-triplet crossover can be observed. While in most cases 
the flows are free of divergences, for parameter sets corresponding to strong triplet 
coupling, the interaction coupling constants diverge for A — >• 0. This divergence 
is probably an artefact of the approximation, in which dynamic and three-particle 
correlations are ignored. Similar divergences appear in some quantum dot geometries 
even in absence of Hund's coupling terms. [2"5] 

Fig.[2]shows the dependence of the conductance for two different values of the level 
splitting, S = 0.5U (curve 1) and 6 = OAU (curve 2), as a function of the gate voltage 
V g which shifts the energy levels according to Eqs. For large 5 one observes two 
well separated conductance plateaus which each arise from the usual spin-1/2 Kondo 
effect. In the first plateau region the occupation of the lowest orbital is one while 
the upper orbital is empty. In the second plateau region the lowest orbital is fully 
occupied whereas the upper orbital is singly occupied. Around V g = the occupancy 
of the lower orbital is close to two (singlet) whereas the upper orbital is almost empty. 
This leads to a strongly suppressed conductance at and around V g — 0. A slightly 
smaller value of the level splitting 5 induces a singlet-triplet transition which leads to 
the occurrence of a conductance plateau with G s=s 4e 2 /h around V g = 0. 

While our approach does not allow to calculate the temperature dependence of 
the conductance, we can investigate its field dependence (at zero temperature) to 
study the stability of the observed conductance plateaus since a small magnetic field 
suppresses the Kondo effect in a smilar manner as a finite temperature. The peak of 
the conductance at V g = is clearly rather fragile, a small magnetic field cuts the 
peak in half (curve 4) while the regular S = 1/2 Kondo plateau hardly changes for 
comparable values of the magnetic field (curve 3). We can in fact define a Kondo 
scale as e.g. the value of Zeeman splitting h produced by a magnetic field necessary 
to suppress conductance by a factor of 1/2. 

If the level separation is made very small and close to V g ~ 0, eventually all 
coupling constants diverge when A approaches zero. The divergence of the parameters, 
which would presumably disappear in a more refined approximation. [23] is however 
such that they define a triplet state. 

As we already mentioned, one of the main advantages of FRG over NRG is the 
ability to reliably calculate zero-frequency properties very quickly. One can thus scan 
large parts of the parameter space to produce 2D maps such as the one in Fig. [3j It 
shows the conductance as a function of V g and Zeeman splitting h. Here, one sees 
again how the triplet plateau is suppressed at a lower value of h than the S — 1/2 
plateau. Since U a = Ub, the conductance is symmetric around V g — and it is thus 
sufficient to plot for V g > 0. 

Figure [4] shows the conduction at the symmetry point V g = as a function of 
the level separation S, for different values of magnetic field H . As seen by the curve 
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Figure 2. Conductance as a function of the gate voltage V g for Ua = Ub — 
U ab = U , J = —0.3(7, Ta = = 0.02ttU. A small change in level separation 
pushes the system into the triplet configuration. (Curves 1 and 2,5 = 0.5U, 0.4(7, 
respectively.) Adding a small magnetic field, one sees that the triplet conduction 
peak is much more fragile than the regular Kondo S = 1/2 plateaus. (1 — > 3, 
2 -> 4.) 



corresponding to H — (curve 1), this regime of 5 covers the singlet-triplet transition. 
Even a very small magnetic field suppresses the conduction on the triplet side strongly, 
but at the transition, a conduction peak remains. The dashed line is an asymptote 
taken from Eq. ( 10 ), valid for large S. The derivation of the asymptotic result Eq. ( 10 1 
does however only take into account the triplet and local singlet configuration of the 
quantum dot and becomes inaccurate for too large 5. As before, we also plot the 
results on a 2D map, shown in Fig. [5j 

To further illustrate the aforementioned triplet divergence of parameters, we plot 
the conductance for V g — as a function of S and J in Fig. [6] The black dots mark 
spots where the flow diverged, clearly showing that this occurs when the system is 
strongly defined to be in the triplet state, with small S and large \J\. The flow of the 
parameters is visualized in Fig. [7] for three sets of initial parameters. The top plot 
shows the flow on the singlet side of the transition, the middle one on the triplet side 
but still close to the transition, and the bottom one far into the triplet side, where the 
divergence occurs. But, as stated before, the parameters diverge in a controlled way: 
| J | and Ua, Ub become very large (J stays negative) while UaB goes to zero, so it is 
clearly preferable to have the two electrons in the triplet state. 



Functional renormalization group approach to the singlet-triplet transition in quantum dotsll 




6. Summary and conclusions 

In summary, we have extended the static FRG approach presented in Ref. [15] for 
multi-orbital quantum dots to include also Hund's coupling and studied the triplet- 
singlet transition in vertically coupled dots, at T — 0. Despite its simplicity, the here 
presented approach reproduces qualitatively the low-temperature results from NRG 
approaches. [5] The dependence of the conductance on the gate voltage follows the 
predictions of the perturbative RG analysis pQ near the triplet-singlet transition. 

Higher order static correlations would probably give rise to small quantitative 
improvements, but a discussion of the full temperature and magnetic field dependence 
of the conductance would require an extension to a dynamic approximation, which, at 
the moment, are only available for the weak to moderate coupling regime. [T2l [TBI H"4] 
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Figure 5. Singlet-triplet transition as a function of h. One can see how the 
conduction on the triplet side gets suppressed and a peak at the singlet-triplet 
transition remains. The swerve of the peak to higher 5 at the top of the plot is 
because the magnetic splitting h is becoming comparable to the x-axis unit. 
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Figure 6. Conductance as a function of <5 and J. Direct Coulomb interactions 
are the standard ones, and we are at the symmetry point V g = 0. As expected 
from Eq. there is a linear relation between <5 and J along the singlet-triplet 
transition curve. However, renormalization effects shift the line such that the 
slope is not 1/4 and it does not go through (5, J) = (0, 0). The slope is about 0.2 
and an extrapolation shows that the curve passes through J = at S ~ 0.37(7. 
The black dots correspond to parameters where the flow diverges. 
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Figure 7. Flow of parameters for Ua = Ub = Uab = U, J = —0.3(7, 
T A = F B = 0.02tt[/, V g = 0. From top to bottom, S = 0.5C7, 0.4(7, 0.38U. 
Top: final value of the inter-orbital Coulomb repulsion is much larger than the 
intra-orbital one, so it the electrons both stay on the bottom level in the singlet 
configuration. Middle: The inter-orbital Coulomb repulsion vanishes and the 
Hund's exchange grows very large, so the triplet configuration becomes favorable. 
Bottom: at very small level separation, the flow diverges at some value of A, but 
at least one can see from the divergence that the system is in the triplet state. 



